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ABSTRACT 

We present an analysis of the 2-point correlation function, £(r), of the X-ray Brightest 
Abell-type Cluster sample (XBACs; Ebeling et al. 1996) and of the cosmological 
constraints that it provides. If £(r) is modelled as a power-law, £(r) = (r /r) 7 , we 
find ro — 26.0 ± 4.5 h~ 1 Mpc and 7 ~ 2.0 ± 0.4, with errors corresponding to 2a 
uncertainties for one significant fitting parameter. As a general feature, £(r) is found 
to remain positive up to r ~ 50-55 /i _1 Mpc , after which it declines and crosses zero. 
Only a marginal increase of the correlation amplitude is found as the flux limit is 
increased from 5 x 10 -12 erg s _1 cm~ 2 to 12 x 10~ 12 erg s _1 cm~ 2 , thus indicating 
a weak dependence of the correlation amplitude on the cluster X-ray luminosity. 
Furthermore, we present a method to predict correlation functions for flux-limited X- 
ray cluster samples from cosmological models. The method is based on the analytical 
recipe by Mo & White (1996) and on an empirical approach to convert cluster fluxes 
into masses. We use a maximum-likelihood method to place constraints on the model 
parameter space from the XBACs £(r). For scale-free primordial spectra, we find 
that the shape parameter of the power spectrum is determined to lie in the 2a range 
0.05^ 0.20. As for the amplitude of the power-spectrum, we find as — 0.4-0.8 for 
Slo = 1 and as — 0.8-2.0 for Slo = 0-3- The latter result is in complete agreement 
with, although less constraining than, results based on the local cluster abundance. 

Key words: cosmology: theory - galaxies: clusters - large-scale structure of Universe. 



1 INTRODUCTION 

Galaxy clusters have been recognised since a long time as 
extremely useful tracers of the large-scale structure of the 
Universe. Being the largest virialized cosmic structures, they 
can be detected rather unambiguously up to very large dis- 
tances. Furthermore, their distribution traces scales which 
have not yet undergone the non-linear phase of gravitational 
clustering, therefore simplifying their connection to the ini- 
tial conditions of cosmic structure formation. 

So far, most of the studies of the cluster distribution 
have relied on optical samples like the Abell/ ACO one (Abell 
1958; Abell Corwin & Olowin 1989), the APM sample (Dal- 
ton et al. 1994) and the EDCC sample (Collins et al. 1995). 
A general result from such analyses is that the two-point 
cluster correlation function is well approximated by a power- 
law, £(r) = (ro/r) 7 , with 7 ~ 1.8-2. The measured correla- 
tion length ro, range in the interval ro — 15-25 /i _1 Mpc , 
depending on the analysed sample and/or on the cluster 



richness (see, e.g., Bahcall & West 1992, Croft et al. 1997; 
and references therein). 

A serious problem arising in optical cluster compilations 
is the spurious enhancement of the cluster richness in high 
density environments which can also cause inherently poor 
clusters or groups to appear rich enough to be included in 
the sample (e.g. van Haarlem, Frenk & White 1997). This 
led several authors to cast doubts about the reliability of the 
clustering analysis of Abell/ACO clusters, whose correlation 
function could be artificially enhanced by such projection ef- 
fects (cf. Sutherland 1988; Dekel et al. 1989; Peacock & West 
1992; see, however, Jing, Plionis & Valdarnini 1992 for an 
alternative view). Although samples like the APM and the 
EDCC have been designed with the purpose of minimising 
projection effects, they could still be present at some level 
(e.g., Collins et al. 1995). 

This calls for the need to confirm results based on opti- 
cal samples by studying other cluster compilations which are 
free of the above mentioned biases. X-ray selected cluster 
samples are ideally suited for this purpose. Indeed, since the 



© 1998 RAS 



2 S.Borgani, M. Plionis, V. Kolokotronis 



X-ray emissivity of the intra-cluster medium (ICM) is pro- 
portional to the square of the gas density, cluster emission 
is strongly peaked at the centre, so as to make the impact of 
projection effects almost negligible. Important analyses of 
the large-scale distribution of relatively small X-ray clus- 
ter samples have already been performed in recent years (cf. 
Romer et al. 1994; Nichol, Briel & Henry 1994, Guzzo et al. 
1995). 

More recently, the X-Ray Brightest Abell Cluster sam- 
ple (XBACs, hereafter; Ebeling et al. 1996) has been com- 
piled from cross-correlating the ROSAT all-sky X-ray sur- 
vey with the Abell/ACO cluster sample (for further details 
see their Section 4). Although it is not purely X-ray se- 
lected, the inclusion of only those Abell/ACO clusters that 
have an extended and significant X-ray emission suppress 
the biases affecting the purely optical sample (cf. Ebeling 
et al. 1996, Plionis & Kolokotronis 1998; PK98, hereafter). 
The XBACs provides for the first time a whole sky, flux- 
limited sample of X-ray galaxy clusters. This sample has 
been analysed by PK98 with the aim of investigating the lo- 
cal acceleration field and by Abadi, Lambas & Muriel (1998) 
to study the cluster 2-point correlation function (see also 
Kolokotronis 1998). 

In this paper, we present a new analysis of the correla- 
tion function for the XBACs cluster distribution, with the 
specific aim of investigating its implication on cosmological 
models for large-scale structure formation. In this context, a 
great advantage in using galaxy clusters, instead of galaxies, 
as tracers of the large-scale structure lies in their less am- 
biguous connection with the underlying dark matter (DM) 
density fluctuations. As for X-ray clusters, the problem of 
understanding their biasing with respect to the DM distri- 
bution requires a suitable description of the connection be- 
tween the X-ray luminosity and the mass of the virialized 
DM halo hosting the emitting gas. This problem, although 
not of trivial solution, is much more affordable, both from 
a numerical and an analytical point of view, than an ac- 
curate understanding of galaxy formation and evolutionary 
processes. 

In what follows, we will use the analytical approach 
by Mo & White (1996) to describe the correlation func- 
tion for the distribution of clusters, which are identified as 
virialized halos of a given mass. This recipe in itself is not 
enough to compare model predictions with XBACs data, 
since the XBACs cluster selection is clearly not based on a 
mass threshold criterion. For this reason, we will resort to 
a phenomelogical approach to connect cluster masses to ob- 
served X-ray fluxes (cf. also Borgani et al. 1998). Although 
our analysis is entirely based on a sample of local clusters, 
all the formalism that we will present is developed to treat 
the cluster distribution at a generic redshift. Therefore, our 
analysis method is amenable to be directly applied to obtain 
cosmological constraints from future deep surveys of galaxy 
clusters. 

The structure of the paper can be summarised as fol- 
lows. In Section 2 we provide a short description of the 
XBACs sample and of the results of the £(r) analysis. We 
describe in Section 3 how to compute cluster correlations 
for an X-ray flux-limited sample from cosmological mod- 
els, which are characterised by a given power-spectrum of 
density fluctuations and by the geometry of the Friedmann 
background. The method to compare model predictions to 



real data analysis is discussed in Section 4. In Section 5 we 
briefly discuss such results and draw the main conclusions. 

2 CORRELATION ANALYSIS OF XBACS 
2.1 The XBACs sample 

The XBACs sample consists of the X-ray brightest 
Abell/ACO clusters that have been detected in the ROSAT 
all sky survey (RASS; Triimper 1990; Voges 1992) for fluxes 
above S Um = 5 x 10" 12 erg s _1 cm -2 ([0.1-2.4] keV band) 
with redshifts limited by z < 0.2. The sample contains 253 
clusters and thus it is the largest X-ray flux-limited cluster 
sample to date (though not entirely X-ray selected). The 
X-ray fluxes measured initially by the Standard Analysis 
Software System (SASS) point source detection algorithm 
are superseded by Voronoi Tessellation Percolation (VTP) 
measurements, which better account for the extended na- 
ture of the emission. In addition, the difficulty of the SASS 
algorithm to actually detect nearby X-ray emission has been 
mostly corrected by running VTP on the RASS fields centred 
on the optical positions of all nearby Abell/ACO clusters 
(z < 0.05) irrespective of whether or not they are detected 
by SASS. Ebeling et al. (1996) estimate, after a careful anal- 
ysis of possible selection effects and biases that the overall 
completeness of this X-ray sample is ^ 80% at the above 
flux limit. 

An extended discussion of the possible XBACs system- 
atic effects is presented in PK98. Here we present only a 
brief summary. 

The known volume incompleteness effect, as a function 
of distance, for the richness class R=0 Abell/ACO clus- 
ters, is not present in the XBACs sample (Ebeling et al. 
1996). Indeed, thanks to its flux-limited property, the sam- 
ple contains at large distances the inherently brighter and 
thus richer Abell/ACO clusters, for which there is no volume 
incompleteness. Furthermore the significant distance depen- 
dent density variations, which are most probably due to the 
higher sensitivity of the ACO Illa-J emulsion plates, be- 
tween the northern Abell and southern ACO parts of the 
combined optical cluster sample (cf. Batuski et al. 1989; 
Scaramella et al. 1991; Plionis & Valdarnini 1991) has been 
corrected for in the XBAC sample (see details in PK98). 

The main remaining biases that need to be quantified 
for the correlation analysis are the redshift selection func- 
tion, which depends on the XBACs luminosity function and 
flux limit, and the Galactic absorption latitude selection 
function, which has been found to be consistent with a cose- 
cant law of the form 

P(\b\) =dex [a(l-csc|6|)] (1) 

with a ~ 0.3 for the Abell sample (Bahcall & Soneira 1983; 
Postman et al. 1989) and a ~ 0.2 for the ACO sample (Ba- 
tuski et al. 1989). In the following we will restrict our anal- 
ysis to clusters having |6| > 30°. 

The XBACs redshift distribution, in the form of an 
N(r) plot, corrected for Galactic absorption, can be seen in 
Figure |jj We have also plotted as a continuous line the pre- 
dicted radial distribution function based on the XBACs lu- 
minosity function (see references in PK98) and the 5 x 10~ 12 
erg s _1 cm -2 flux limit. Note that redshifts are converted 
into luminosity distances according to 
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Figure 1. The XBACs radial distribution function and its Pois- 
son uncertainty (filled circles). The predicted distribution based 
on the XBACs luminosity function estimated by PK98 is shown 
as a continuous line. 
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where Hq — 100 h km s— 1 Mpc 1 is the Hubble constant 
and 



'(*) = 



dz E 1 (z) ; Qa = 1 — fio 



2 \n z + (2 - fio) (1 - VI + itoz)] 

Here £(2) = [(1 + 2) 3 fi + n A ] 1/2 for a flat Universe (e.g., 
Peebles 1993) . In what follows, results and comparisons with 
models will be presented for two values of the density pa- 
rameter; fio = 1 and 0.3, with and without a cosmological 
constant term (fiA = A/3-Hq) to restore spatial flatness. 



2.2 The £(r) estimate 

We estimate the 2-point correlation function for XBACs us- 
ing the following estimator: 



l + £(r) =2 



(N cr ) 



(4) 



where N cc is the number of cluster pairs in the interval 
[r — 8r,r + Sr] and (N cr ) is the average, over 10,000 random 
simulations with the same boundaries, redshift and galac- 
tic latitude selection functions, cluster-random pairs in the 
same separation interval. We have evaluated f (r) in logarith- 
mic intervals with constant logarithmic amplitude A ~ 0.1. 
We estimate the variance of £(r) by using the analytical 
approximation to the bootstrap errors (Mo, Jing & Borner 
1992) which has been shown to reproduce fairly accurately 
the actual bootstrap variance: 



Table 1. Results for the three samples. Column 2: flux limits 
(units of 10 — 12 erg s — 1 cm -2 ); Column 3: number of clusters; Col- 
umn 4: correlation length and 2a uncertainty (units of h~ 1 Mpc ); 
Column 5: slope of £(r) and 2a uncertainty. 



(a) 



No. of clusters 



203 



(b) 8 

(c) 12 



112 



(i7 



ro 



26.CT 



25.9 -1 



1.98 



+0.35 



27.5 



+8.7 



2.02 



1.94"* 



+0.43 



3.9 



3.6 



5.8 



2 



3 x 



2 



= 3 X 



N cc (r) 



(5) 



where a 2 P is the quasi-Poisson variance. 

We apply the correlation analysis to the whole sam- 



ple (z< 0.2), with Siim 



5 x 10" 12 erg s 1 cm" 2 



sample 



(a)] and to two other subsamples with Sn m = 8 x 10" erg 
s" 1 cm" 2 and SW = 12 x 10" 12 erg s" 1 cm" 2 [samples (b) 
and (c), respectively]. The total number of clusters in each 
sample is reported in Table [l]. 

The resulting £(»"), using Qo — 1 in e< p([|), for the three 
samples are plotted in Figure ^. Only for sample (a) we 
show the £(r) by using fio = 0.3. Due to the limited depth 
of XBACs no significant differences are found between these 
two cases. For all the three samples £(r) is positive up to 
r ~ 50 h, _1 Mpc , after which it declines and crosses zero at 
about 50 - 55 /i -1 Mpc . This result agrees with that found 
for optical Abell/ACO cluster samples (e.g., Klypin & Rhee 
1994) and confirms that the cluster distribution requires a 
substantial amount of large-scale coherence of cosmic den- 
sity fluctuations. 

The dashed lines in the plots represent the best-fitting 
power law model, £(r) = (ro/r) 7 , which is determined 
through a x 2_m inimisation (Xmin hereafter) procedure, as- 
suming a Gaussian distribution for o^. The fit for the (c) 
sample has been performed only including bins with r > 10 
h~ Mpc , so as to exclude the anomalous high-correlation 
signal coming from the smallest scales. On a more quantita- 
tive ground, Figure ^ shows the iso-A^ 2 contours (A% 2 = 
X 2 ~ Xmin> with Xmin being the absolute minimum value of 
X 2 ) in the ro — f plane. Strictly speaking, £(r) estimates in 
different bins are not independent of each other, since each 
cluster contributes to pairs at different separations. The con- 
tours correspond to lcr and 2er uncertainties for two signif- 
icant parameters and correspond to A^ 2 = 2.30 and 6.17, 
respectively. Figure ^] shows the variation of A^ 2 around the 
best-fitting value of each of the two parameters, once we 
marginalise with respect to the other parameter. The best 
fitting values for ro and 7, along with the 2cr uncertainty for 
one significant fitting parameter are also reported in Table 
|l| It turns out that increasing Sum the correlation length 
only marginally increases. As long as a correlation exists 
between cluster A^-ray luminosity and richness, this result 
implies only a mild dependence of the clustering strength on 
the cluster richness (see also Croft et al. 1997). 

Results for sample (a) can be compared with those ob- 
tained by Abadi et al. (1998) for the same sample. They 
find a somewhat lower correlation amplitude, with ro — 21 
/i -1 Mpc and a similar 7 value but with much smaller uncer- 
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Figure 2. The 2-point correlation function for samples (a), (b) and (c), by assuming f2o = 1 to convert redshifts into distances (cf. 
eqs.(^) and (^). Errorbars correspond to 1<t bootstrap uncertainties. Open circles for the (a) sample corresponds to the open Qo = 0.3 
case. The dashed lines corresponds to the best fitting power-law expressions, §(r) = (ro/r) 7 , whose parameters are reported in Table ^. 



tainties (cf. their Fig. 4) . We are not sure why their estimate 
of ro is lower than ours. However, the rather large difference 
in the uncertainties is entirely due to their assumption of 
Poissonian errors for £(r). We verified that repeating our 
analysis assuming quasi-Poisson errors (cf. the upper right 
panel in Fig. ^), the contours in the ro~7 plane shrinks into 
a much smaller size, comparable to that reported by Abadi 
et al. (1998). 



3 PREDICTING CORRELATIONS FOR 
FLUX-LIMITED SAMPLES 

In this section we introduce the formalism to compute model 
£(r) for a flux-limited cluster survey. After briefly introduc- 
ing the analytical method by Mo & White (1996) to estimate 
correlations for virialized halos of a given mass, we present 
an empirical procedure to convert this mass into an X-ray 
flux in the appropriate energy band. A similar approach has 
been also applied by Moscardini et al. (1999, in prepara- 
tion) to predict the 2-point correlation functions expected 
in different X-ray flux-limited cluster samples. We point 
out that, although our method is applied here to the analy- 
sis of the local XBACs clusters, it is presented in a general 
form so as to be directly applicable to any sample of both lo- 
cal and distant X-ray selected clusters. Therefore, it is also 
well suited to study the clustering evolution using future 
deep cluster surveys. 

3.1 The analytical recipe for cluster correlations 

Our estimates of model cluster correlations is based on the 
approach originally proposed by Mo & White (1996; MW 
hereafter) . In this approach, the correlation function at red- 
shift 2 for clusters (identified as virialized halos) of mass 



M, £ c i(r, z, M), is connected to the dark matter correlation 
function at the same redshift, £_dm(»", z), according to 



where 

(,DM(r,z) = 



,M)^dm{t,z) 



\D{x) 




[D(0)_ 





dkP(k) k 



sin kr 
kr 



(6) 



(7) 



In the above expression D(z) is the linear growth factor for 
density fluctuations at redshift z (e.g., Peebles 1993). 

The biasing factor, b(M, z), appearing in eq.(^) is given 

by 



KM, 



1 + 



[8 c {z)/g{Rm, 



S c (z) 



(8) 



(e.g, MW; Matarrese et al. 1997; Baugh et al. 1998). Here 
S c (z) is the critical linear overdensity for spherical collapse. 
For a critical-density Universe 8 C = 1.686, independent of z, 
with a weak dependence of the geometry of the Friedmann 
background (e.g., Eke, Cole & Frenk 1996). Furthermore, 



a{R 



M,Z) = 



\D( Z y 




[D(0)_ 


2^ 2 Jo 



(9) 



is the fluctuation variance at mass M and redshift z, with 
W(x) the Fourier representation of the window function, 
which describes the shape of the volume from where the col- 
lapsing object accreates matter. The comoving fluctuation 
size Rm represents the Lagrangian cluster radius. It is con- 
nected to the mass scale M and to the present day average 
matter density p according to 



Ra 



( Mi 

I 47rp 



1/3 



(10) 



© 1998 RAS, MNRAS 000, |-?? 



Cosmological Constraints from the £(r) of the XBACs 5 



i 1 1 1 i 



S??=5 




I I I I I I I II 



i i i i i i i i 



50 
40 
30 
20 
llO 

a 

£.50 
^40 
30 
20 
10 



Figure 3. Iso— Ax 2 contours on the rg-j parameter space from 
the correlation analysis of samples (a), (b) and (c) (upper left, 
lower left and lower right panels, respectively). The upper right 
panel show the results for sample (a) if quasi-Poissonian errors 
for f (r) were instead assumed. Contours correspond to la, 2a and 
3a confidence levels for 2 significant fitting parameters. 
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parameter. Continuous, dotted and dashed curves are for (a), (b) 

and (c) samples, respectively. The two horizontal lines and the 

upper limit of each panel represent the la, 2a and 3a confidence 

levels for one significant parameter. 



for the top-hat window, W(x) = 3(sinx — :rcos:r)/:r 3 , that 
we adopt in the following. 

Since in practical cases one is interested in the correla- 
tion of halos with mass above a given limit M, the biasing 
factor in eq.(^) should be replaced by an effective bias, & e //i 
whose value is obtained by averaging b over the mass distri- 
bution of the virialized halos: 



be ff (M,z) = 



dM'b(M',z)n(M', 
Jm dM' n(M',z) 



In the above expression the mass distribution n(M, z) is 
given by the Press & Schechter (1974) expression 



n(M,z)dM = 
S 2 c (z) 



exp 



2a 2 (M,z) 



o 8 c (z) 

7T M 2 a (M, Z ) 

dM . 



dloga(M, z) 



dlosM 



(12) 



The power-spectrum P(k) of density fluctuations can be 
expressed as P(k) = Ak n T 2 (k), where n is the primordial 
spectral index. A Harrison-Zel'dovich primordial spectrum, 
with n — 1, will be assumed in the following. As for the 
transfer function, T(k), we take the expression 



= md + 2.34,) 
v; 2.34<? 

[l + 3.89q + (16.1g) 2 + (5.46q) 3 + (6.71g) 4 ] 



-1/4 



(13) 



(11) 



Here q — k/hV, Y being the shape parameter. For the class 
of CDM models, it is F ~ n h (Bardeen et al. 1986), while 
in general V can be viewed as a free parameter, to be fitted 
against observational constraints. For instance, the power- 
spectrum of APM galaxies provide F = 0.23 ± 0.04 (e.g., 
Peacock & Dodds 1996). As for the amplitude of P(k), it is 
customary to express it in terms of as, the r.m.s. fluctuation 
amplitude within a top-hat sphere of 8 /i -1 Mpc radius. 

Therefore, each model for large-scale structure for- 
mation will be characterised by three parameters, namely 
ilo, o"8 and F. In the following, we will present results on the 
r-<78 plane, fixing the density parameter to either Qq = 0.3 
or Qo = 1- 

The reliability of the MW approach to describe the clus- 
tering of cluster-sized halos has been already tested by Mo, 
Jing & White (1996; cf. also Governato et al. 1998). Jing 
(1998) recently showed that a correction to eq.(^|) is required 
only for halos of mass much smaller than M* (defined as 
the mass for which Sc/cm = 1), a regime which is not rel- 
evant for clusters in plausible cosmological models. Colberg 
et al. (1998) found from the analysis of their Hubble Volume 
Simulations that the MW method overpredicts the cluster 
correlation length ro by ~ 20%. In the following analysis 
we will compute the biasing factor according to eq.(|ll]) by 
assuming for 5 C the canonical value of the spherical top-hat 
collapse model. We just point out that, since 6 e // depends 
on S c and cm only through their ratio, a change in the choice 
of 5c turns into a proportional change in the resulting values 
of <t 8 ■ 

Finally, since the real data are analysed in redshift 
space, we introduce the effect of z-distortion in the MW 
expression for Redshift-space correlations are ampli- 

fied by the usual factor K(J3) = 1 + 2/3/3 + (3 2 /5 (Kaiser 
1987), where (3 ~ 0§" 6 /6 e //- However, for the cluster case 
where /3 ~ 0.20, the effect is rather small since it increases 
£(r) by only ~ 15%. 



3.2 Converting fluxes into masses 

Cluster masses are the input quantities for the MW analyt- 
ical approach that we have just described. Since the observ- 
able quantity is the cluster luminosity rather than its mass 
we should provide a suitable method to convert mass into 
fluxes (or vice- versa). 
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As a first step we convert masses into X-ray tempera- 
tures. According to the spherical collapse model and under 
the assumption of virial equilibrium, the mass-temperature 
relation can be written as 



k B T 



1.38 

Pt 



M 



2/3 



10 15 /i-iM© 
n A„ ir (z)] 1/3 (1 + z)keV , 



(14) 



where 76% of the gas is assumed to be hydrogen (see, e.g., 
Eke et al. 1996). A V i r (z) is given by the ratio between the 
mean density within the virialized region and the average 
cosmic density at redshift z. For Slo = 1 it is Aj, ir = 187r 2 , 
independent of z, while we use the expressions provided by 
Kitayama & Suto (1996) for the Qo = 0.3 cases. For isother- 
mal gas, the J3t parameter is defined as the ratio of the spe- 
cific kinetic energy of the collisionless matter to the specific 

2 

thermal energy of the gas, /3r = ^ffi* > with /i = 0.59 the 
mean molecular weight, m p the proton mass and a v the one- 
dimensional cluster internal velocity dispersion. Equivalence 
between specific gas thermal energy and DM kinetic energy 
implies /3t = 1. However, neither of such assumptions may 
be completely correct. 

The calibration of the fir value using numerical simula- 
tions has been attempted by several authors (see, e.g., Bryan 
& Norman 1998, for a summary of numerical results). Re- 
cent simulations of the Santa Barbara Cluster Comparison 
Project (Frenk et al. 1998), based on a variety of numeri- 
cal techniques, indicate that pr — 1.15. This value will be 
adopted in the following as the fiducial one, while we will 
also show the effect of changing this parameter over a rather 
broad range. For instance, a smaller /3t gives a smaller mass 
and, therefore, a larger ctm, at a fixed temperature, so as 
to reduce the biasing factor in eq.(|i"l|). One should however 
bear in mind that such N-body calibrations of /3r generally 
include only adiabatic physics of the intra-cluster medium, 
while neglecting the effects of radiative cooling, feedback ef- 
fects from galactic winds, etc. 

As a second step we convert X-ray temperatures into 
bolometric luminosities, Lbol- This is probably the most del- 
icate step of our procedure, since neither observations nor 
theoretical modelling converge to a unique, well determined 
Lboi-T relation. By adopting a phenomenological approach, 
we model the Lboi~T relation as 



Lboi = L% 



2 -1 

erg s 



(15) 



where Le gives the expected luminosity for a 6 keV cluster, 
while A defines the redshift evolution of the Lb i~T relation. 
Data for local clusters indicates that L$ ~ 3 as a rather 
stable result, and a ~ 2.5-3.5, at least for temperatures 
1 keV, depending on the sample, data analysis tech- 
nique and corrections for cooling flow effects. White, Jones 
& Forman (1997) analysed a set of 207 EINSTEIN clusters 
and found a ~ 3, thus in agreement with previous results 
(e.g., David et al. 1993, and references therein). Although 
the formal fitting uncertainties are generally small, the scat- 
ter of data points around the relation (jl5|) is so large as to 
raise the question of whether it represents a good model for 
the observational Lb i~T relation. A remarkable reduction 
of the scatter is found once the effects of cooling flow are 
properly introduced, at least for temperatures fesT 1 ^ 4keV 



(e.g., White et al. 1997; Arnaud & Evrard 1997; Markevitch 
1998; Allen & Fabian 1998). In the following we will adopt 
the value a = 3 as the reference one, while we will also show 
the effect of its variation on the resulting model constraints. 
Since available data of the redshift dependence of the Lboi~T 
relation indicate no evolution up to z ~ 0.4 (Mushotzky & 
Scharf 1997) we will use A = in the following. No depen- 
dence of the final results on A are expected, owing to the 
limited depth of XBACs. 

As a third step we convert bolometric luminosities, Ltoi , 
into finite energy-band ([0.1-2.4] keV) luminosities, Le, ac- 
cording to Le = /e[T(M, z); z]Lboi, where the bolometric 
correction term, /b, is computed by integrating the emis- 
sivity over the appropriate energy band. Following Math- 
iesen & Evrard (1998), we assume in our analysis a purely 
Bremsstrahlung emissivity, with the power-law approxima- 
tion for the Gaunt factor, g[E,T) oc {E /kT)~ ' 3 . This ap- 
proximation has been shown to be rather precise for TxzL 2 
keV, which is the temperature range expected to be relevant 
for the bright Abell clusters of XBACs, in correcting lumi- 
nosities in the soft ROSAT band, [0.5-2.0] keV (Borgani et 
al. 1998). We expect this approximation to be at least as 
good, when dealing with the broader XBACs band. 



4 ANALYSIS AND RESULTS 

In order to place constraints on the model parameter space, 
we should device a method which is able to estimate the val- 
ues of the slope F of the power spectrum and its amplitude 
as, along with the respective statistical uncertainties, in the 
most assumption-free way. 

Binning the data as in Fig. |i| may in principle introduce 
an uncertainty in the estimated £(r), whose value depends in 
general on the binning choice. Furthermore, the theoretical 
£(r) of eq.(|^) requires the value of the mass threshold M to 
be known. Even for a fixed choice of the mass-luminosity 
relation, the value of M depends on the effective redshift, 
z e ffi that has to be associated to the whole survey. Indeed, 
for a given flux-limit, the effective luminosity limit is given 
by 



Leff = 4:TYd L (z e ff)Slir, 



(16) 



where (Ll is the luminosity-distance. In principle, a reason- 
able choice for z e ff is thepeak of N(z), which occurs at 
~ 0.075 - 0.08 (see Fig. ^). However, due to the broad- 
ness of the XBACs redshift distribution, the precise value 
of z e ff can not be fixed in an unambiguous way. We show 
in Figure |B| the effect of changing the value of z e ff from 
0.05 to 0.1 on £(r), for a fixed cosmological models and 
Lx-Tx relation. Note that, for both such values of z e ff, 
the redshift distribution is far from declining and, there- 
fore, any redshift within this interval can in principle be 
adopted as z e ff. Increasing z e ff from 0.05 to 0.1, has a non- 
negligible effect on the resulting f(r): the larger z e ff, the 
larger the luminosity- and, therefore, the mass-threshold. In 
turn, this corresponds to a larger biasing factor and, hence, 
to an increase of £(r). Indeed, at z e ff = 0.05, the mass- 
threshold is M t h — 2.1 x W 14 ^ 1 Mq which corresponds to 
an effective biasing 6 e // ~ 4.6. Such numbers increases to 
Mth ^ 4.4 x 10 14 /i _1 M Q and b eff ~ 6.2 for z eff = 0.1. 

This demonstrates that any ambiguity in the choice of 
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Figure 5. The influence of the choice for the effective survey red- 
shift, z e ff, on the cluster 2-point correlation function predicted 
by a specific cosmological model with a fixed Lx-Tx relation. 
The flux— limit corresponds to that of the sample (a), whose, £(r) 
is also plotted as filled circles. 

z e ff generates a similar ambiguity in the determination of 
the best-fitting model. 

Furthermore, at large separations the estimate of £(r) is 
likely to be dominated by more distant and, therefore, more 
luminous clusters. As a consequence, £(r) at different scales 
takes contributions from cluster populations having different 
luminosity and, in principle, different clustering properties. 
A possible solution to this problem would be selecting only 
clusters above a fixed luminosity limit. It is however clear 
that this remedy can only be pursued at the expense of 
substantially reducing the sample statistics. 

To overcome such limitations of the binned Xmin ap- 
proach, one can estimate the likelihood that the model cor- 
relation function produces the measured number of cluster 
pairs at a given separation and at a given redshift, for spec- 
ified sample flux-limit and mass-luminosity relation, in a 
way which does not depend on the binning procedure. 

Let £(r, z, M(z)) be the model correlation function for 
halos that, at redshift z, have mass above the value required 
by the sample flux-limit. Then, the number of cluster- 
cluster pairs with separation between r — dr/2 and r + dr/2 
and at redshift between z — dz/2 and z + dz/2 is 

C{r,z)drdz = [1 +£ c i(r,z,M(z))] 1*^1*1 drdz. (17) 

Here lZ(r, z) is the expected number of cluster-random 
pairs, within the same separation and redshift interval. We 
estimated TZ(r) by averaging over 20,000 random samples. 
The two-dimensional grid in r and z has 50 equal amplitude 
bins in the separation interval 6 < r < 56 ft Mpc and 10 
equal amplitude bins for redshift z < 0.2. It is clear that the 
bin size can be made as small as desired, by proportionally 
increasing the number of random samples, so as to make the 
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Figure 6. Iso-AS contours on the ag-F plane from the sample 
(a), for fio = 1 and 0.3, with both flat and open geometries. 
The upper right panel show the results from the ^-minimisation 
analysis for the critical-density case. Contours correspond to la, 
2a and 3it confidence levels for 2 significant fitting parameters. In 
each panel, the shaded area represents the la confidence region 
from the COBE normalisation, as given by White & Scott (1997). 

final results independent of the bin size. The relevant quan- 
tity to estimate is the likelihood function C, which is denned 
as the product of the probabilities of having one pair at each 
of the (ri, Zi) bins occupied by the data cluster-cluster pairs 
and the probability of having zero pairs in all the other el- 
ements of the r—z plane. By assuming Poisson probabilities 
for bin occupation, we get 

C(r,z) = TT exp[— C(r, z) drdz] C(r, z) dr dz 

i 

x Y[exp[-C{r,z)drdz] , (18) 

where the indices i and j runs over the occupied and empty 
bins, respectively. The best fitting values of the model pa- 
rameters are found by minimising the quantity S = — 21n£. 
After dropping all the terms independent of the cosmological 
scenario, it reads 

S = 2 J dr J dzC(r,z)-2^2lnC(n,Zi) . (19) 

i 

A similar maximum-likelihood (ML hereafter) approach has 
been applied by Croft et al. (1997) to estimate the richness 
dependence of the correlation length ro of APM clusters (see 
also Marshall et al. 1983). 

The results of the application of the ML analysis are 
shown in Figure ^, for a low-density model with S7o = 0.3, 
both with open and fiat geometry, and for the critical den- 
sity case. While the selected range of the shape parameter 
r turns out to be almost independent of f2o, the value of as 
is a decreasing function of the density parameter, much like 
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the as values coming from the local cluster abundance con- 
straint (e.g., Eke et al. 1996; Viana & Liddle 1996; Girardi 
et al. 1998). It is worth remembering that the confidence 
regions in the as~T plane from the ML analysis are esti- 
mated with the underlying assumption that cluster pairs are 
independent of each other. A careful error estimate would 
pass through the construction of a non-Gaussian likelihood 
function (e.g., Dodelson, Hui & Jaffe 1997), which in princi- 
ple requires the knowledge of the whole TV-point correlation 
function hierarchy. In their ML analysis of APM clusters, 
Croft et al. (1998) compared errors from the Poissonian like- 
lihood with the cosmic variance from N-body simulations. 
They concluded that accuracy of errors from a Poissonian 
ML function depends on the cluster richness, being rather 
accurate for cluster populations having a mean separation 
d c i^ 70 /i -1 Mpc . Assuming a — 3 in the Lx~Tx relation 
and Pt = 1.15 in the M-Tx conversion, the XBACs flux- 
limit corresponds to a cluster mass-limit ~ 3.8 x 10 14 fe _1 M© 
at z e ff — 0.08. In turn, taking the Qo = 1 best-fitting model 
with (78 = 0.55 and F = 0.11, the resulting cluster abun- 
dance from the Press-Schechter formula would correspond 
to d c i ~ 85 fe _1 Mpc for the mean separation of clusters 
above this mass limit. Therefore, we expect the confidence 
regions found by the ML approach to provide a realistic esti- 
mate of the actual uncertainties in the parameter estimates. 

In order to check the robustness of the ML results, we 
also decided to apply the Xmin analysis to the binned £(r). 
As already mentioned, a potential limitation of the Xmin 
approach lies in the ambiguity in the definition of the effec- 
tive redshift of the survey. In this analysis, we decided to 
assume z e ff = 0.075, which corresponds to the peak of the 
redshift distribution. The results for the Xmin analysis for 
Oo = 1 are reported in the upper right right panel of Figure 
^[ It turns out that the two methods identify essentially the 
same regions of the as~T plane. From the one hand, this re- 
sult supports the robustness of the analysis methods. From 
the other hand, it indicates that the effective survey redshift 
is in this case reliably represented by the peak of N(z). We 
report in Table 2 the best-fitting parameter values, for the 
three cosmological models, for both the ML and the Xmin 
method, along with their 2a uncertainties for one significant 
parameter, i.e. after maginalising with respect to the other 
parameter. 

We note that as takes larger values for the smaller 
Qq. This is due to the fact that, for fixed cluster distance 
and flux-mass conversion, the sample flux limit turns into a 
larger value of the cluster Lagrangian radius R m in eq.(|lo|) 
when a smaller density parameter is assumed. Therefore, for 
a fixed P(k) amplitude, clusters would correspond to more 
rare and thus more clustered peaks. A suitable increase of as 
(i.e., of the P(k) amplitude) is thus required to decrease the 
relative height of the peaks to be associated with clusters. 

The resulting constraint on as and its dependence on 
Qo is comfortably consistent with, although somewhat less 
stringent than, that coming from the abundance of local 
galaxy clusters. Indeed, several independent analyses based 
on the optical virial mass function (e.g., Girardi et al. 1998), 
the local X-ray temperature function (e.g., Viana & Liddle 
1996; Eke et al. 1996; Oukbir, Bartlett & Blanchard 1997) 
and the local X-ray luminosity function (e.g., Borgani et 
al. 1998) converge to indicate that as — 0.5-0.6 for Slo = 1, 
while it rises to as — 0.9-1.2 for Oo = 0.3, almost inde- 




Figure 7. The variation of AS around the best-fitting value of 
each parameter, after marginalising with respect to the other pa- 
rameter. Continuous, dotted and dashed curves are for the (a), 
(b) and (c) samples, respectively. As in Figure ^, the two hori- 
zontal lines and the upper limit of each panel represent the la, 
2a and 3a confidence levels for one significant parameter. 
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Figure 8. The same as Figure |jj but for the effect of changing 
the parameter a in the Lx—Tx relation. Dotted, continuous and 
dashed curves are for a = 2.5, 3 and 3.5, respectively. 

pendent of the presence of a cosmological constant term. 
This shows that combining results from cluster abundance 
and clustering do not allow to break the degeneracy be- 
tween as and Q.q. In any case, it is remarkable that both 
the abundance and the large-scale distribution of clusters, 
which involve rather different scale ranges, can be consis- 
tently interpreted in the framework of hierarchical clustering 
of Gaussian density fluctuations, as described by the Press 
& Schechter (1974) approach and by its MW extension to 
the 2-point correlation function. 

As for the shape parameter, it is generally constrained 
to lie in the range 0.05,$ F$ 0.20, thus in agreement with, al- 
though somewhat smaller than, that obtained from the opti- 
cal Abell/ACO cluster distribution (Borgani et al. 1997) and 
from the APM galaxy distribution (e.g., Peacock & Dodds 
1994). 

The shaded areas in each panel of Fig. |^ represent the 
la constraint from the COBE normalisation, as provided 
by White & Bunn (1996). It turns out that, if fi = 1, 
as — 0.4-0.5 and F ~ 0.18 (at la c.l.) are required to satisfy 
at the same time the large-scale CMB constraint and the 
XBACs cluster clustering. As for the low-density flat model, 
it requires as — 0.9-1.1, with essentially the same shape 
parameter. As for the open model, CMB and cluster £(r) 
constraints are rather inconsistent with each other. 

It is worth reminding here that the above results refer 
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Figure 9. The same as Figure [?], but for the effect of changing 
the mass-temperature conversion factor, /3t- Dotted, continuous 
and dashed curves are for 0t = 0.8, 1.15 and 1.5, respectively. 



to the whole XBACs, with Su m = 5 x 10~ 12 erg s" 
well as to one particular choice of the Lx~Tx relation and 
of the M-Tx conversion factor Pt- 

In Figure (?) we plot for the fio = 1 case the variation 
of AS = S — S m in around its minimum, as one of the two 
parameters as and F is varied, keeping the other fixed at 
its best fitting value, for the (a), (b) and (c) samples. It 
turns out that results are quite stable when varying the flux 
limit of the XBACs sample analysed; the only significant 
effect being that of enlarging the confidence intervals as the 
number of selected clusters decreases with increasing Su m . 
The stability of the results against variations of the flux limit 
indicates that any possible incompleteness of XBACs at low 
fluxes has at most a marginal effect. 

The effect of changing the shape of the local Lx~Tx 
relation is shown in Figure pi The extreme values of a = 2.5 
and 3.5 have been chosen so as to bracket the range of results 
presented in the literature (cf. the discussion in the previous 
section). Even in this case, constraints on model parameters 
are rather robust against changes in the Lx~Tx relation. 
This represents a remarkable result, since connecting lumi- 
nosities to temperatures represents in principle a rather del- 
icate step, owing to the uncertainties in both observational 
and theoretical determinations of the Lx~Tx relation. 

Finally, we show in Figure |^ the effect of changing the 
parameter /3t in the M-Tx relation. Even in this case, the 
extreme values, /3t = 0.8 and 1.5 have been chosen so as 
to largely encompass the results from hydrodynamic cluster 
simulations (cf. Bryan & Norman 1998). We find that, while 
the shape parameter F does not depend on /3t, the value 
of as systematically increases with /3t- The reason for this 
is that a larger /3t implies a larger mass for a fixed clus- 
ter temperature. As a consequence, a larger as is required 
to suppress the biasing factor, so as to compensate for the 
selection of more massive, more clustered clusters. 



5 CONCLUSIONS 

We presented an analysis of the redshift -space 2-point cor- 
relation function for the X-ray Brightest Abell-type Cluster 
sample (XBACs; Ebeling et al. 1996), with the aim of inves- 
tigating the resulting constraints on cosmological models. 
In our analysis, cosmological models are specified by three 
parameters, namely the density parameter fio (with or with- 



Table 2. Model parameters for the sample (a) from the maximum 
likelihood and from the x 2 -uiinimisation analyses. Uncertainties 
corresponds to 2er confidence levels for one significant parameter. 



Max-like. 



X — min. 



Uo 

n 



i 

0.3; U A 
0.3; U A 



CT8 

0.55 





0.7 



1.07 
1.12 



+0.23 
0.11 
+0.88 
0.27 
+0.88 
0.29 



0.11 
0.11 
0.13 



+0.05 
0.03 
+0.06 
0.04 
+0.06 
0.05 



CT8 

0.59 



1.28 
1.50 



0.30 
0.12 
+0.96 
0.26 
+ 1.33 
0.47 



0.11 

0.11 
0.11 



0.07 

0.04 
+0.12 

0.06 
+0.13 

0.06 



out a cosmological constant term to provide spatial flatness), 
the fluctuation power-spectrum amplitude through the as 
quantity and the shape F of the power-spectrum. 

As a starting point, we follow the analytical approach 
by Mo & White (1996), which provides the 2-point corre- 
lation function, £(r), for the distribution of virialized halos 
above a given mass limit (e.g., Mo, Jing & White 1996). In 
order to convert masses into observed fluxes, we followed a 
purely empirical recipe (see also Borgani et al. 1998), which 
depends on two parameters, namely the conversion factor 
(3t to pass from mass to X-ray temperatures [cf. eq.(n4j)] 
and the slope a of the local Lx~Tx relation [cf. eq.([15j)~ 

Although this analysis is limited to nearby (z^ 0.2) 
clusters, the formalism that we introduced can be directly 
applied to any survey of high-redshift objects. 

The comparison between model and data has been per- 
formed by resorting both to a x 2- niinimisation (Xmin) pro- 
cedure and to a maximum likelihood (ML) approach, which 
avoids any ambiguity associated to the binning procedure 
required by the Xmin method. 

The main results of our analysis can be summarised as 
follows. 

(a) If the 2-point correlation function is modelled as a 
power-law, £(r) = (ro/r) 7 , then the best fitting parame- 
ters for the whole sample are ro = 26.0 ± 4.5 ft _1 Mpc and 
7 = 2.0 ± 0.4. The clustering strength increases very weakly 
as a higher flux-limit is imposed and higher luminosity clus- 
ters are selected. As long as a correlation exists between clus- 
ter luminosity and richness, this result is consistent with the 
picture of a mild dependence of the correlation amplitude on 
the cluster richness (cf. Croft et al. 1997). 

(b) For our reference choice for the mass-luminosity con- 
version (i.e., a = 3 and @t = 1-15), we find that as — 
0.6 ± 0.2 (2a uncertainties after marginalisation) , while it 
becomes about twice as large for £Iq = 0.3. This result is 
completely consistent with constraints coming from cluster 
abundances and indicates that the picture of hierarchical 
clustering of Gaussian primordial density fluctuations is able 
to account at the same time for both the abundance and the 
clustering of galaxy clusters. 

As for the shape of P(k) it is much less dependent on flo, 
with values ranging in the interval 0.05,$ T^S 0.20. Such re- 
sults are left unchanged by increasing the sample flux-limit. 

(c) Adding also the large-scale CMB constraints, we find 
that an no = 1 model with as — 0.4-0.5 and F ~ 0.16- 
0.20 and a flat fio = 0.3 model with as — 0.9-1.1 and the 
same F are consistent with both the COBE data and the 
clustering of XBACs clusters. For the open flo = 0.3 model, 
the COBE and XBACs constraints are quite inconsistent 
with each other. 

(d) As for the robustness of such results against changes 
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of the mass-luminosity relation, it turns out that F is quite 
insensitive to both the mass-temperature conversion factor 
(3t and to the slope a of the local Lx-Tx relation. The 
power spectrum amplitude erg has only a weak dependence 
on the choice for the Lx~Tx relation, while it increases from 
about 0.35 to 0.6 as /3t varies from 0.8 to 1.5. 

As a concluding remark, we point out that, although 
XBACs provides useful constraints on cosmological mod- 
els, it i a not a completely X ray aclcctcd aamplc. Newer 



and lar gor X - ray cluetor eamploe havo boon rocontly com — 

piled (BCS; Ebeling et al. 1998) or are close to completion 
(REFLEX; Bohringer et al. 1995). Such samples, thanks to 
their careful definition of completeness criteria and to their 
increased statistics will allow to definitely fix the cluster- 
ing scenario for local X-ray clusters. Even more exciting, 
the possibility of realizing cluster surveys at higher redshifts 
with next-generation X-ray satellites will provide a further 
means to constrain both the power spectrum of density fluc- 
tuations and the value of the matter density parameter. 
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